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In Hanbury-Brown-Twiss interferometry measurements using identical bosons, the chaoticity pa- 
rameter A has been introduced phenomenologically to represent the momentum correlation function 
at zero relative momentum. It is useful to study an exactly solvable problem in which the A param- 
eter and its dependence on the coherence properties of the boson system can be worked out in great 
detail. We are therefore motivated to study the state of a gas of noninteracting identical bosons at 
various temperatures held together in a harmonic oscillator potential that arises either externally 
or from bosons' own mean fields. We determine the degree of Bose-Einstein condensation and its 
momentum correlation function as a function of the attributes of the boson environment. The pa- 
rameter A can then be evaluated from the momentum correlation function. We find that the X(p, T) 
parameter is a sensitive function of both the average pair momentum p and the temperature T, and 
the occurrence of A = 1 is not a consistent measure of the absence of a coherent condensate fraction. 
In particular, for large values of p, the A parameter attains the value of unity even for significantly 
coherent systems with large condensate fractions. We find that if a pion system maintains a static 
equilibrium within its mean field, and if it contains a root-mean-squared radius, a pion number, 
and a temperature typical of those in high-energy heavy-ion collisions, then it will contain a large 
fraction of the Bose-Einstein pion condensate. 

PACS numbers: 25.75.-q 25.75.Gz 03.75.-b 05.30.Jp 

I. INTRODUCTION 

In high-energy collision processes, Hanbury-Brown-Twiss (HBT) intensity interferometry [1] has been used to study 
the space-time geometry of the source of particles [2] -[38]. As is well known, for identical bosons the interference 
phenomenon arises from Bose-Einstein correlations and depends sensitively on the degree of coherence of the boson 
system [39]. The HBT correlation occurs for a chaotic source but not for a coherent source. 

In phenomcnological measurements, one represents the correlation in terms of the momentum correlation function 
C(p + q/2,p-q/2) = C(p 1 ,p 2 ) = G (2 \p l ,p 2 ;p l ,p 2 )/G ( - 1 \p ll p l )G^{p 2 ,p 2 ), where p x and p 2 are the momenta of 
the pion pair and G^' is the i-th-order pion density matrix. One introduces the parameter A = [C(q = 0) — 1] that is 
purported to represent the degree of chaoticity of the pion medium and bears the name "the chaoticity parameter" . 
Experimental measurements with pions persistently indicate that the chaoticity parameter A is substantially less than 
the value of unity for a fully chaotic source. Some part of this reduction of the chaoticity parameter A from unity 
may be attributed to the occurrence of the decays of long- live resonances [40] . However as emphasized not the least 
by Glauber [41], part of the reduction of the 'chaoticity parameter' A from unity may arise from the coherence of the 
pion gas. 

Even though the chaoticity parameter has been widely used in all HBT measurements in high-energy collisions, how 
the chaoticity parameter can be determined theoretically has not been resolved. The difficulty is further compounded 
for heavy-ion collisions because the dynamics of pions after their production in high-energy heavy-ion collisions is 
very complicated. The process of initial hadronization and the subsequent interactions between pions are beyond the 
realm of present-day knowledge. 

It is therefore useful at this stage to study an exactly solvable problem for which the A parameter can be determined 
explicitly and the transition from the coherent phase to the chaotic phase can be worked out in detail. We are 
motivated to investigate the state of a gas of noninteracting identical bosons held together in a harmonic oscillator 
potential at various temperatures. We shall study the occurrence of Bose-Einstein condensation and the two-body 
momentum correlation function as a function of the attributes of the boson gas in such an environment. This will 
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allow us to examine explicitly the transition from the coherent phase to the chaotic phase and to study how this phase 
transition may affect the HBT measurements and the A parameter for a set of known attributes of the Bose-Einstein 
gas assembly. 

In atomic physics, the harmonic oscillator potential introduced here can arise from an external trap. In high- 
energy heavy-ion collisions, the harmonic oscillator potential can arise approximately from the mean-field potential 
experienced by a pion, owing to the interactions generated by other pions and medium particles. Although the 
strength of the pion mean field is not known at present, the results obtained here will serve as useful supplementary 
tools to study the circumstances in which a pion system may form a Bose-Einstein condensate in heavy- ion collisions. 
They will stimulate future investigations on the magnitude of the pion mean-field potential and pave the way for 
future investigations on momentum correlations for pions under more complicated dynamical evolutions. 

Pions produced in high-energy heavy-ion collisions have a temperature that is of the order of the pion rest mass. 
The motion of the pions is relativistic and a proper treatment will need to be relativistic in nature. We shall carry 
out both a nonrelativistic and a relativistic treatment of the pions to understand what features of the coherence are 
sensitively affected by the relativistic motion. 

Important advances in our understanding of the coherence properties of identical bosons have been made recently 
in another related field, the physics of atomic boson systems at low temperatures [42-51]. Theoretical and experi- 
mental work in atomic physics has focused on the correlation function in the configuration space. In particular, the 
second-order correlation function (r±, r^) has been obtained to give the probability of detecting a boson at r\ 
in coincidence with the detection of another identical boson at r 2 [43] . From the shape of this correlation function 
gr( 2 )(n, r%) as a function of the relative separation r x — r 2 , the theoretical spatial correlation length can be extracted. 
Experimentally, the measurements of various arrival times and positions at the detectors in HBT intcrfcromctry are 
then used to determine the spatial correlation length, for comparison with theoretical analyses [50] . We wish to adopt 
a treatment complementary to that in atomic physics by examining the correlation function in momentum space, the 
standard arena for Bose-Einstein correlation analysis in high-energy nuclear collisions [2]- [38]. Our investigation of 
the correlation function in momentum space is greatly facilitated by utilizing the results of the correlation function 
in configuration space obtained in atomic physics [42-44] . 

With regard to the low-temperature measurements with atoms, the perspectives of studying the correlation in 
momentum space presented here offer useful complementary points of view. In momentum space the trapped atoms 
are now described as having an equilibrium momentum distribution, appropriate for the system in a given external 
field at a given temperature. The sudden removal the external field allows the initial momentum distribution of 
the particle to be frozen at the moment of the external field removal, as appropriate under the application of the 
sudden approximation in quantum mechanics. Subsequent free streaming of the particles without the external field 
and mutual interactions allows the reconstruction of the momentum distribution of the source at the moment of 
its freezing out. In atomic physics, the correlation function in the complementary momentum space has many rich 
features as it is sensitive to many kinematic variables and the geometry of the source particles. 

This paper is organized as follows. In Section II, we review the degree of Bose-Einstein condensation as a function of 
temperature and particle number. In Section III, we study the one-body and two-body momentum density matrices. 
In Section IV, we express the momentum correlation function in terms of the one-body momentum density matrix and 
the ground-state wave function and also express the momentum correlation function in terms of the Wigner function. 
In Section V, we evaluate the one-body density matrix and the Wigner function for bosons in the harmonic oscillator 
potential. In Section VI, we study the spatial and momentum distributions of these boson assemblies. In Section VII, 
we evaluate the momentum correlation function C(p,q), for different values of the average pair momentum p and 
temperature T. We extract the A parameter and the HBT radii. In Section VIII we study the condensate fraction for 
a nonrelativistic pion gas with a given root-mean-squared radius in a static equilibrium at various temperatures. In 
Section IX, we investigate the relativistic treatment of the boson in a harmonic oscillator potential. In Section X, we 
evaluate the boson spatial density and estimate its condensate fraction for the relativistic boson gas. In Section XI, 
we present our conclusions and discussions. 

II. CONDENSATE FRACTION AS A FUNCTION OF TEMPERATURE 

We shall first review the theoretical work on the coherence and correlations of identical bosons in atomic physics 
[42-44] so as to pave the way for our investigation of the correlation function in momentum space. We consider first 
a nonrelativistic gas of identical bosons in a harmonic potential at temperature T with the potential specified by 



in 



3 



where m is the rest mass of a boson and hv measures the st rength of the external potential. We shall measure lengths 
in units of the harmonic oscillator length parameter a = \JhJmuj, momenta in units of h/a, and energies in units of 
huj. 

The states in the harmonic oscillator potential are characterized by energy levels e„ = (n+ |)?zu; with the associated 
degeneracy of g n = (n+ l)(n + 2)/2. Following Ref. [43], it is convenient to use the recalibrated energy level e„ = nha; 
measured relative to 3H/2. 

As the temperature of the gas is lowered below the condensation temperature T c , condensation of the nonrelativistic 
massive boson gas occurs. As is well known, the fluctuation of the number of particles in the condensate state, the 
n = state, depends sensitively on the assumed statistical ensemble. A grand canonical ensemble will lead to a 
condensate ground state number fluctuation that is as large as the number of particles in the condensate ground 
state, n — 0. The grand canonical ensemble cannot be used to describe the number of particles Nq in the ground 
state condensate. The condensation can best be studied in a canonical ensemble for the case with a fixed number of 
particles [42]. Comparison of the results from the canonical and the grand canonical ensemble in Ref. [42] indicates 
however that even though the number of particles in the ground n — state can only be described by the canonical 
ensemble, the occupation number distribution of the n > harmonic oscillator states can be appropriately described 
by the grand canonical ensemble with only very small corrections. The difference in these two ensembles for the n > 
states becomes very small as the number of particles increases. 

Therefore, for a fixed number of particles iV at a given temperature T/Jko, we shall follow Rcfs. [42] and [43] to 
determine the condensate configuration by the following requirements: (i) a fixed total number of particles TV in a 
canonical ensemble for the condensate n = state, and (ii) an occupation number distribution in a grand canonical 
ensemble for the n > states. Accordingly, we have the following three conditions to determine the condensate 
configuration of the system with a fixed number of particles N at a temperature T = 1/(3: 

N = N + N T , (2) 

where N is the number of condensate particles in the n = state, 

N ° = (3) 
Nt is the number of "chaotic" particles in the n > states, 
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and z is the fugacity parameter. For the harmonic oscillator potential, the summation for Nt can be carried out 
analytically and Eq. (4) can be simplified to 

, T fe e" fe ^(3- s e -k^ + e -2k0^ 

Nt = ^ Z (i_ e -^)3 ( 5 ) 

Equations (2)-(4) can be reduced into a single condensate configuration condition, 

z ^ k e~ kl3hul (3 - 3 e -k/3hw + e -2k/3hujj 
fe=i 



iV-^+Vz fc ^ K - — '- (6) 



Because N and f3huj = fku/T are fixed, the above condensate configuration condition can be solved numerically to 
determine the unknown z (by Newton's method with fast convergence). After the value of the solution z is obtained, 
iVo and Nt can be subsequently determined from Eqs. (3) and (5) to give the condensate configuration specified by 
the condensate fraction f and the 'chaotic' fraction fx, 

N N T 
fo = — and h = w . (7) 

We show in Fig. 1 the fugacity solution z which satisfies the condensate configuration condition Eq. (6) for different 
temperatures T/fou) and boson numbers N. To get a better view of the z values, we show an expanded view of Fig. 
1(a) in the z <~ 1 region in Fig. 1(6). We observe that the fugacity parameter z is close to unity in the strongly 
coherent region at low temperatures. In fact, the fugacity parameter z at T = assumes the value 
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FIG. 1: (Color online) (a) The fugacity parameter z satisfying the condensate configuration condition Eq. (6) for different 
boson numbers N, as a function of temperature T/hui and (6) an expanded view in the z ~ 1 region. 

For a given boson number N, as the temperature increases from T — 0, the fugacity z decreases very slowly in the 
form of a plateau until the condensate temperature T c is reached, and it decreases very rapidly thereafter. The greater 
the number of bosons N, the greater is the plateau region, as shown in Fig. 1(6). For example, for N = 2000 the value 
of z is close to unity for < T/hu; < 11 in the plateau, and it deviates from unity substantially only for temperatures 
T/tiw» 11. 




FIG. 2: (Color online) Solid curves represent the condensate fractions fa(T), calculated with the condensate configuration 
condition Eq. (6), as a function of T/hio for different boson numbers N. The abscissa labels for the corresponding chaotic 
fraction fr(T) = [1 — fo(T)] are indicated on the right. The dashed curves are the fits to the solid curve results of fo(T) with 
the function 1 — (T/T c ) 3 of Eq. (10) where the values of T c /hu> for different N values are listed in Table I. 

The condensate fractions fo(T) calculated with the fugacity parameter of Fig. 1 for different boson numbers N 
are represented by solid curves in Fig. 2, as a function of T/tko. The abscissa labels for the corresponding chaotic 
fraction fr(T) = [1 — fo{T)] are indicated on the right. We observe that the condensate fractions are unity at T = 0, 
corresponding to a completely coherent boson system at T = 0. It decreases slowly as the temperature increases, and 
the rate of decrease is small at low temperatures. The greater the number of bosons N, the larger is the range of 
temperatures in which the boson system contains a substantial fraction of the condensate. For example, for a system 
with 2000 identical bosons, substantial fraction of the condensate occurs up to T/hw ~ 11. The transition from the 
condensate phase to the chaotic phase occurs over a large range of temperatures and is therefore not a sharp first-order 
type transition. The complementary chaotic fraction fr(T) increases gradually as the temperature increases, reaching 
the value of unity at T/Hu ~ 11 for N = 2000. 

As the number of particles N decreases down to 250, a substantial condensate fraction occurs only for T/hu < 5. 
The chaotic fraction Jt(T) increases as a function of temperature and it reaches the value of unity for T/fko ~ 6. 
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The transition from the condensate phase to the chaotic phase occurs over a temperature range from T/Hco ~ 2 to 
T/ftw ~ 5. The smaller the number of particles, the lower the condensate temperature T c and the smaller is the range 
of temperatures over which the condensate phase transition occurs. 

In the transitional region below T c with a substantial fraction of the condensate, one can get an approximate value 
of the condensate fraction by noting that in this region, the value of z is close to unity (Fig. 1). The number of 
chaotic particles Nt can be estimated from Eq. (5) by setting z to unity, and we obtain 

(rji \ 3 OO 
7 fc=i 

Consequently one can fit the condensate fraction fo(T) reasonably well by a one-parameter function of the form 

f (T) = 1 - (T/T c f for T<T C , (10) 



1 2{hw/T) 15(fej/T)' 



(9) 



f (T) = 0(1/N) - for T>T C 



(11) 



The results from the one-parameter fit to fo(T) are shown as the dashed curves in Fig. 2, to be compared with the 
/o(T) calculated with the condensate configuration condition Eq. (6) shown as the solid curves. The values of T c /hw 
that give the best fit to fo{T) for different N values are listed in Table I. 

The T c values can also be determined approximately by considering the case of fko/T « 1 in Eq. (9) and we have 
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where £(3) = X^fcLi ^ 3 = 1-202 is the zeta function with the argument 3. Thus, the condensate fraction is given 
approximately by 



f (T) ~ 1 - (T/T c . approx f for T<T C , 



with 



c.approx 
fojj 



( N 



VC(3) 



1/3 



N 
1.202 



1/3 



(13) 



(14) 



A comparison of the above approximate result with T c in Table I indicates that the above Eqs. (13) and (14) are 
approximately valid, with the values of T c approx determined by Eq. (14) slightly greater than T c by about 10%. 

Table I. Condensation temperature T c /hui of Eq. (10) 



and T, 



c.approx 



/Huj of Eq. (14) as a function of N 



Number of Bosons N 


T c /tlLU 


T c , approx /hw = (A7l.202) 1 /3 


^c.approx i T c 


2000 


10.97 


11.85 


1.08 


1000 


8.56 


9.41 


1.10 


500 


6.63 


7.47 


1.13 


250 


5.12 


5.92 


1.16 



III. ONE-BODY AND TWO-BODY DENSITY MATRICES IN MOMENTUM SPACE 

Previously, the one- and two-body density matrices have been obtained for identical bosons in configuration space 
[43]. We would like to write down the corresponding one-body and two-body density matrices in momentum space 
so as to evaluate the momentum correlation function. The results in momentum space can be readily obtained from 
the results in configuration space by replacing r in Ref. [43] with p. We thus have the one-body density matrix in 
momentum space 

G (1) (Pi,P 2 ) =E !1 n(Pl) M °(P2)(«i«n). (15) 
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Similarly, we have the two-body density matrix in momentum space given by 

G (2) (Pl,P 2 ;Pl,P2) = Wk(PlX(P2)Wm(p 2 )Wn(Pl)(4«l t ama n ). (16) 

klmn 

We shall follow [43] in expressing the two-body density matrix in terms of one-body density matrices. By separating 
out the term with k = 1 = m = n from other terms and using the definition of the one-body density matrix (15), the 
two-body density matrix can be shown to be 

G (2) ( Pl ,p 2 ; Pl ,p 2 ) - GW(p 1 ,p 1 )G( 1 )(p 2 ,p 2 ) + |G«(p 1 ,p 2 )| 2 

oc 

+ ^|<( Pl )| 2 |u n (p 2 )| 2 {(at4a n a n )-2(ata n )(ata n )}. (17) 

n=0 

The last term in this equation involves a summation over the n = condensate state and the set of {n > 0} states. 
In line with our earlier discussions on the statistical ensemble for the states [42, 43], we shall use the grand canonical 
ensemble for the set of {n > 0} states and the canonical ensemble for the condensate state of n = 0. For the 
set of {n > 0} states in the grand canonical ensemble then, the occupation fluctuation characteristics of the grand 
canonical ensemble make the contributions of the set of {n > 0} states small in comparison with the other terms on 
the right-hand side of the above equation, as we shall see from the following discussion. We note that in the above 
equation 

(44^^) - 2(ala n )(al l a n ) = ((a£a n - (a^a n )) 2 ) - (aj 1 d„)(a] 1 a n ). (18) 

For an n > state in the grand canonical ensemble, the mean-square fluctuation of the occupation number (a£a n ) in 
the state is given by [52] 

((4 a n - (ata n )) 2 ) = (ala n ) ((a£a n ) + l) . (19) 
Therefore, for this n > state in the grand canonical ensemble, we have 

(alala n a n ) - 2(o£a n )(a| 1 a„) = (ala n } 7 (20) 
and the contribution of the set of {n > 0} states to the two-body density matrix is 

^|<( Pl )| 2 | Un (p 2 )| 2 (at an ). (21) 

n>0 

When we integrate over p l and p 2 , the set of {n > 0} states gives a contribution of 

d Pl dp 2 K(Pi)| 2 K(p 2 )| 2 <a£a n > - N T , (22) 

n>0 

whereas the other terms such as the first term of Eq. (17), ( Pl , p 1 )G^ 1 - ) (p 2 , p 2 ), give a contribution of TV 2 . The 
contribution from the set of {n > 0} states is Nt/N 2 of the contribution from G^ 1 \p 1 ,p 1 )G ( - 1 \p 2 ,p 2 ). Therefore, 
in the limit of a large number of bosons N, the ratio Nt/N 2 is small, and the contributions from the set of {n > 0} 
states in the summation in Eq. (17) can be neglected. We are left with only the n = condensate state contribution 
for this summation. 

To describe the n = condensate state, we shall follow Ref. [42, 43] and use the canonical ensemble which gives 
the canonical fluctuation [42] 

((a£a n - {ala n )) 2 ) = (ajajaoao) - (aJa )(aJao) = O(N ). (23) 

Thus, we have 

{alala n a n ) - 2(ala n )(al a n ) = -(aJa )(aJao) + O(N ). (24) 

In the limit of a large number of particles, we can neglect the last term O(iVb) in the above equation which is small 
compared to the first term of order Nq. The two-body momentum density matrix of Eq. (17) is therefore 

G (2) ( Pl ,p 2 ; Pl ,p 2 ) - G^( Pl , Pl )G^(p 2 ,p 2 ) + \G (1 \ Pl ,p 2 )\ 2 - N 2 \u ( Pl )\ 2 \u (p 2 )\ 2 , (25) 

which gives the conditional probability for the occurrence of a pion of momentum p 1 in coincidence with another 
identical pion of momentum p 2 . This two-body density matrix in momentum space has the same form as that 
obtained earlier in configuration space in [43]. 
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IV. THE MOMENTUM CORRELATION FUNCTION 



In Hanbury-Brown-Twiss measurements, we normalize the probability relative to the probability of detecting particle 
p l and p 2l and define the momentum correlation function C(p 1 ,p 2 ) as 

c( v u)- g(2) (pi^ 2 ;pi 1 p 2 ) (26] 
C(Pl ' P2) -Ga)( Pl , Pl )Ga)(p 2) p 2 )- (26) 

It is convenient to introduce the average and the relative momenta of the pair 

P=(Pi+P 2 )A q = Pi~P 2 , (27) 

with the inverse transformation 

Pi=P+2' P2=P-7}- (28) 

The momentum correlation function can be expressed alternatively in terms of the kinematic variables p and q. From 
Eq. (25), we have the general expression for the correlation function 

In the nearly completely coherent case with almost all particles in the ground condensate state, No — > N, the two 
terms in the numerator cancel each other and we have C(p,q) = 1, as it should be. For the other extreme of a 
completely chaotic source with No << N, the second term in the numerator proportional to Nq gives negligible 
contribution and can be neglected. The correlation function becomes the usual one for a completely chaotic source, 

r r„^ |g (1) (Pi,P2)l 2 

C chaotic (p, q) = 1 + mi —^ m -— y (30) 

The general result of Eq. (29) allows one to study the correlation function for all cases with varying degrees of 
coherence. 

If we introduce R{p 1 q) = R(p 1 ,p 1 ) = C(p, q) — 1, then 

Mr, a) -Mr, n) ^Wgjj 2 - ^o 2 ko(Pl)| 2 ko(P 2 )| 2 

K{p,q) - K{p 1 ,p 2 ) — rp- \ ■ ( 61 > 

G {1 >(Pl,Pl)G {1, (P2,P2) 

It is of interest to express the momentum correlation function C(p, q) in terms of the Wigner function f(r,p) defined 
in terms of the the one-body density matrix G^(ri, r 2 ) as 

f(r,p) = J dsjP s G^{r+ S -,r-^), (32) 

This one-body density matrix in configurations space is related to the one-body density matrix in momentum space 
by a Fourier transform, 

G (1) (Pi,P 2 ) = / dri dr 2 e^-^G^ (r,, r 2 ). (33) 

Therefore, by changing coordinates from n and r 2 to r — (n + r 2 )/2 and s = n — r 2l we can relate the one-body 
density G' 1 ^(p 1 ,p 2 ) with the Wigner function f(r,p), 

G (1) (Pi,P 2 )- J dre^ r f(r,p), (34) 

and in particular, for the diagonal density matrix element we have 

G (1) (Pi,Pi) = J drf(r, Pl ). (35) 
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As a consequence, the correlation function C(p, q) in Eq. (29) can be rewritten as 

r( \ i . | J dre^ r /(r,p)| 2 - Nj\uo(p + g/2)|> (p - g/2)| 2 

{P ' q) ~ + fdrf(r,p + q/2)fdrf(r,p-q/2) ' W 

This is the general expression for the momentum correlation function expressed in terms of the Wigner function 
f(r,p) when the coherence of the system is properly taken into account. 

The R function for the general case is related to the Wigner function f(r,p) by, 

r>r n) = \Jdre^ r /(r,p)| 2 - N$\u (p + g/2)| 2 | Mo (p - g/2)| 2 
(P,q> fdrf(r,p + q/2)fdrf(r,p~q/2) ' {6 '> 

When the condensate fraction f is large with iV — > N, the second term in the numerator of the above equations 
is important and must be properly taken into account. In fact, in the completely coherent case, Eqs. (36) and (37) 
give C(p, q) = 1 and R(p, q) = 0. Only in the special case of a completely chaotic state is the contribution from the 
second term in the numerator negligible, and we have the usual relationship between the Wigner function and the 
momentum correlation function for a chaotic system, 

I f dre i( l r f(r oil 2 

C chaotic (p, q) ~ 1 + f ,'-V /9 , /; 'P ToV ( 38 ) 

J drf(r,p + q/2) J drf(r,p-q/2) 



V. THE ONE-BODY DENSITY MATRIX FOR A HARMONIC OSCILLATOR POTENTIAL 

For a given total particle number N of particle mass m in an external harmonic oscillator potential, we have 
obtained in Section II the fugacity z as a function of T/fou) (Fig. 1). This solution of z allows us to evaluate the 
density matrices and the correlation functions at various temperatures. For the harmonic oscillator potential, the 
one-body density matrix has been obtained previously in configuration space as given by [43] 

ze -0e a 



G«(n,r 2 ) = £<(riK(r 2 ) _ 

n=0 

oo 

= J2 zk Go(ri,r 2 ;kf3fru), (39) 
k=l 

where 

r, , x ( 1 \ 3/2 ( 1 (r 2 + r 2 )(cosh T-l) + (n-r 2 ) 2 \ 
G ° (ri ' r2;r ^ U(l- e - 2 -) j CX H'^ 2sinhr (40) 

Because of the exchange symmetry of r/a and pa/h for a harmonic oscillator potential, the one-body density matrix 
in momentum space can be readily obtained from these results of Ref. [43] by replacing r/a with pa/h, and we get 

ze -0e» 



n=0 Ze 

OO 

= 5> fe G ( Pl ,p 2 ;fc/3M, (41) 



k=l 

where 



Go(Pl ' P2 ' T) - U^l-e-^J P 2^inh7 ) • 

We can write Gq(p 1 ,p 2 ;t) m terms of the ground state wave function u^(p 1 )u (p 2 ) as 



Go(Pi,P 2 ;t) = u* (p 1 )u (p 2 )g (p 1 ,p 2 ;T) 



(43) 
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where the ground state wave function is 



rf>= J r l h> 2 



and the dimensionless function go{Pi,P 2 \ T ) is given by 

1 / a 2 (p\ +p^)(cosh t - 1 - sinhr) + (p 1 -p 2 ) 2 \ 
^P^) = (T^T^exp [~ ¥ y -^^ 2sinhr j" (45) 

Then we have 

G (1) (Pi,P 2 ) = K(Pi)MP2)MPi,P2) (46) 

where 

oo 

A{ Pl ,p 2 ) = J2 zk 9o(Pi,Pi;kf3hiu). (47) 

fe=i 

In numerical calculations, especially at low temperatures where z is close to unity, the number of terms in the 
summation over k in A(p 1 ,p 2 ) will need to be greater than the number of particles N. To avoid such a lengthy 
summation, it is simplest to separate out the condensate component to write the above as 

oo 

A( Pl , p 2 ) = JL- +^z k [g {p 1 , p 1 ; kphu) - 1] . (48) 
z k=i 

For low temperatures, the coefficient [go(Pi,Pi', kf3huj) — 1] of z k is small and a small number of terms in k will suffice. 
For high temperatures above the condensate temperature, z is substantially less than unity, and z k decreases rapidly 
as k increases; a small number of terms in k will also suffice. 
From Eq. (29) the momentum correlation function is 

and 

T>( \ T3( n ) |A( P1 , P2 )| 2 -|Z/(1~Z)| 2 

The one-body Wigner function for the boson system can be obtained from the one-body density matrix G^\ri,r 2 ) 
and we find 

/<-» -±* mmT^{- (s + *£) - m\ 

There is an explicit symmetry between x/a and pa in the Wigner function for the harmonic oscillator potential. 



VI. SPATIAL AND MOMENTUM DISTRIBUTIONS 



Before we evaluate the momentum correlation function C (p, q) , it is useful to study the single-particle spatial and 
momentum distributions p r (r) and p p (p). Because of the symmetry between r/a and pa/h in a harmonic oscillator 
potential, the following two functions have the same shape: (i) p r (r) in units of a~ 3 expressed as a function of r/a, 
and (ii) p p (p) in units of (a/H) 3 expressed as a function of pa/h. The two distributions can be displayed on the same 
graph. From the one-body density matrix (46), we obtain 

P P (p) = GW(p,p)= ^^ 3/2 CX p|-^|yl(p,p), (52) 
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r/a , path 



FIG. 3: (Color online) The spatial density distribution p r {r) in units of a -3 , expressed as a function of r/a, and the momentum 
density distribution p P (p) in units of (a/ft) 3 , expressed as a function of pa/h, for a boson system with TV = 2000 at different 
temperatures. 

where 

oo 
k=l 

1 f a 2 p 2 (cosh t — 1 — sinh r) \ 
5o(P,P;r) = (1 _ e _ 2T)3/2 exp (-- — j . (54) 

We plot in Fig. 3 the spatial and momentum distributions of the system with N — 2000 as a function of their 
dimensionless variables r/a and pa/h, respectively. One observes that up to T/hw ~ 9 the system has a small 
spatial or momentum size and there is a substantial condensate fraction in the system. In Fig. 4 we plot the root- 
mean-squared radius in units of a, r rms /a = \J ((r/a) 2 ), and the root-mean-squared momentum in units of h/a, 
Prmsi/fi- = v ({po^/K) 2 ), as a function of T/hu>. For TV = 2000, the quantity r rms /a is slightly greater than 1 up to 
T/Hlo ~ 6, and it increases relatively rapidly to about 5.5 at the condensate temperature, T c /fuo ~ 11. It increases at 
a relatively slower rate at temperatures above T c . 

The size of the momentum distribution also undergoes similar changes as a function of temperature. The root- 
mean-squared momentum has the dimension of about one unit of h/a at T/hu> ~ and this linear size increases about 
sixfold when the temperature reaches the chaotic region of T c /hui ~ 11 for N = 2000. 

We observe therefore that for a boson system in a harmonic oscillator, the Bose-Einstein condensation gives rise to 
a distribution localized in the region of small momentum and small spatial coordinates. From the viewpoints of the 
spatial and momentum densities, the Bose-Einstein condensate in a harmonic oscillator is therefore a "condensation" 
in both momentum space and configuration space. 

VII. EVALUATION OF THE MOMENTUM CORRELATION FUNCTION C(p, q) 

With the solution z obtained for a given T/hto as shown in Fig. 1 and discussed in Section II, one can use Eq. 
(48) to evaluate A(p l ,p 2 ). The knowledge of A(p l ,p 2 ) then allows the determination of the momentum correlation 
function from C(p,q) using Eq. (49). 

We show results of C(p,q) for the case of N = 2000 in Fig. 5. We observe that the correlation function is a 
complicated function of the average pair momentum p and temperature T. For p = h/a in Fig. 5(a), the correlation 
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FIG. 4: (Color online) The root-mean-squared radius in unit of a and the root-mean-squared momentum in units of h/a, as a 
function of T/hu for different numbers of bosons in the system. 
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FIG. 5: (Color online) The correlation function at different values of the pair momentum pa/h and temperatures. Panels (a), 
5(b), and (c) are for p = 1, 2, and 3 h/a, respectively. 



function C(p,q) at q — is close to unity for temperatures below and up to T/fkv = 9, but increases to 2 rather 
abruptly at T/huj — 12. For p — 2h/a in Fig. 5(6), the correlation function C(p,q) at q = is substantially above 
unity and increases gradually as temperature increases. For p = 3H/a in Fig. 5(c), the correlation function C(p, q) at 
q = is about 2 for all cases of temperatures examined. 

If one follows the standard phenomenological analysis and introduces the chaoticity parameter A to represent the 
correlation function at zero relative momentum, then this parameter A is a function of the average pair momentum p 
and temperature T 

\(p,T) = [C(p,q = 0;T)-l], (55) 

where we display explicitly the dependence of the correlation function on the temperature T. We plot the values 
of X(p,T) as a function of p in Fig. 6(a) for different temperatures for the case of N — 2000. At T/hui = 12, 
which is above the condensate temperature T c , the A parameter is 1 for all p values. At T/fiio = 9, which is below the 
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FIG. 6: (Color Online) (o) The parameter A as a function of pa/h for different temperatures for N 
Rhbt/o, as a function of pa/h for different temperatures for N = 2000. 



2000. (b) The ratio 



condensate temperature T c , the A parameter drops precipitously to ~0.1 at pa/h = 1. At this T/Hlo — 9, as p increases 
the A parameter rises gradually and reaches the constant value of 1 at pa/h — 2.4. At T/Tilu = 6 and 3, for which the 
systems are significantly coherent with large condensate fractions, the A parameter starts close to zero at pa/h — 1, 
but as p increases the A parameter increases gradually to unity at pa/h = 2.9 and 3.1 for T/hio = 6 and 3 respectively. 
The location where the A parameter attains unity changes with temperature. The lower the temperature, the greater 
is the value of p at which the A parameter attains unity. 

We conclude from our results that the parameter A(p, T) is a sensitive function of both p and T and X(p, T) = 1 is 
not a consistent measure of the absence of the condensate fraction, as it attains the value of unity in some kinematic 
regions for significantly coherent systems with large condensate fractions at temperatures much below T c . Only for 
the region of small p will the parameter X(p,T) be correlated with, but not equal to, the chaotic fraction Jt(T) of 
the system. 




FIG. 7: (Color online) (a) The ratio Rhbt/o, as a function of pa/h for different temperatures for iV = 250. (b) The parameter 
A as a function of pa/h for different temperatures for TV = 250. 
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One can evaluate the root-mean-squared momentum in the correlation function defined as 

a 2 (y T) - (a 2 ) - f^BP^L^ m 

q llas {p,l)-{q ) - J dg[c{pq . T) _ l} (5b) 
and introduce the HBT radius -Rhbt (p, T) defined by 

The HBT radius -Rhbt(p, T) is in fact the radius parameter in the standard Gaussian parametrization of the momen- 
tum correlation function, 

C(p, q;T) = 1 + \(p, T) exp{- g 2 i? 2 BT (p, T)/h 2 }. (58) 

We plot Rhbt(p,T) as a function of pa/h for different temperatures T for the case of N — 2000 in Fig. 6(6). 
For fixed values of temperature T/hw — 3,6, and 9 and varying p, one observes that Rbbt(p,T) is about 1.3a 
and Rhbt(p,T) decreases slightly before it increases gradually as p increases. For T/hw = 12, which is above the 
condensate temperature, -RhbtCp, T) is about 2.4a at pa/h = 1, and Rb:bt(p,T) decreases slightly before it increases 
gradually as p increases. For a fixed value of p at pa/h = 1, the HBT radius increases with increasing temperatures 
very slowly at low temperatures, and it increases rather abruptly when the condensation temperature is approached. 
For this small value of pa/h, the variation of Rbbt(p, T) as a function of T reflects closely the variation of the root- 
mean-squared radius as a function of T/hu, as shown in Fig. 4. For large values of p, Rbbt{p,T) increases with 
increasing temperatures in a more uniform manner. 

In Fig. 7(a) and 7(6), we show, respectively, \{p,T) and Rhbt{p,T) as a function of pa/h and T for the case with 
N = 250. The A parameter and the HBT radius -Rhbt behave in a manner similar to those for the case of N = 2000. 
One observes in Fig. 7(a) that at temperatures below T c , the A parameter is small at small p and it increases as p 
increases, reaching the saturating value of unity at pa/h — 3.5 for T/hw = 2 and at pa/h — 2.7 for T/hw = 4. Above 
the condensation temperature at T/hw ~ 6, the A parameter assumes the value of unity for all p values. 

As shown in Fig. 7(6), for a fixed value of temperature T, the HBT radius -Rhbt decreases slightly and then 
increases gradually with p. For a fixed p with a small p, the increase in -Rhbt is slow at low temperatures and the 
increase becomes more rapid as the temperature approaches the condensate temperature of T c /hw = 5.12. 



VIII. BOSE-EINSTEIN CONDENSATION OF PIONS IN A MEAN FIELD (NONRELATIVISTIC) 

There is not much information on the magnitude of the mean-field potential experienced by the pions. From the 
Glauber theory [53], the mean-field potential experienced by a pion in a pion medium is related to the pion density 
p r (r) by 

v{T)= _»m {T) (59) 

m 

where /(0) is the forward ir-ir scattering amplitude. We hope to evaluate the pion mean-field potential in the future. 
In the meantime, the results in the previous sections allow us to answer the following theoretical question. If a 
system of N pions is held together by its mean field, taken to be a harmonic oscillator, and if it comes to a state of 
static equilibrium with a given root-mean-squared radius r rms at a temperature T, what is the condensate fraction 
of such a system? The answer to this theoretical question will provide useful information on the importance of the 
Bose-Einstein condensation for a pion system in static equilibrium, to pave the way for future investigations for the 
system in dynamical expansion. 

We would like to examine pion systems with a typical r rms , T and the pion number that one encounters in high- 
energy heavy-ion collisions. For a pion gas distribution with an HBT radius of about -R HBT = 6 fm as appropriate 
for Au-Au central collisions [9], the root- mean-squared radius r rms for a Gaussian density distribution is v^-Rhbt, 
which is about 10 fm. We shall therefore examine a pion system with r rms = 10 fm, a temperature range from 80 
to 160 MeV, and the number of pions to be N — 250 (for a central SPS Au-Au collision at ^/snn — 19.4 GeV) and 
N = 2000 (for a central RHIC Au-Au collisions at ^/IW^ = 200 GeV). 

For a pion with a temperature of 80 to 160 MeV which is of the order of the pion rest mass of 140 MeV, the 
motion of the pions is relativistic and the proper treatment will need to be relativistic in nature. We shall carry out 
a relativistic treatment of the pion states in the next section and shall content ourselves here in the type of solution 
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one gets in a nonrelativistic treatment. Carrying out both relativistic and nonrelativistic treatments will allow one to 
understand what features of the coherence are sensitively affected by the relativistic motion. 

We first determine the strength of the mean-field potential ftw that can hold a system of N pions in static equilibrium 
at temperature T for a given root-mean-squared radius of r rms . For the pion system in static equilibrium, the quantity 
r rms /a is a function Fn(x) of the variable x = T/huj as shown in Fig. 4, where the subscript N labels the boson 
number. If the value of r rms is fixed as given, the quantities fko and T are then related by the set of parametric 
equations 

^ = l^g, (60) 



T = x 



[hF N (x)f 



(61) 



By varying x = T/ftw for a fixed r rms and using the function F^{x) of Fig. 4 in the above equations, the energy Huj 
can be determined as a function of T for the cases of N = 2000 and TV = 250. The results are shown in Fig. 8(a). 
One finds that for the pion system with a given root-mean-squared radius of 10 fm, the value of fuv ranges from about 
12 to 20 MeV for N = 2000 and about 20 to 30 MeV for N = 250. The ratio of T/fuv is about 7 for N = 2000, 
and is about 4.5 for N = 250, as shown in Fig. 8(6). From these ratios of T/Huj, one can use Fig. 2 to find out 
the condensate fraction. The condensate fractions fo(T) for a pion gas at various temperatures with N = 2000 and 
N = 250 are shown in Fig. 8(c). One finds that f {T) is about 0.67-0.8 for N = 2000 and is about 0.9 for N = 250. 
The knowledge of hui in Fig. 8(a) allows one to determine the values of a as a function of the temperature as shown 
in Fig. 8(d). The oscillator length a is about 4 fm for N = 2000 and about 3.5 fm for N = 250. 
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FIG. 8: (Color online) (a) the potential strength hui, (b) the ratio T/hui, (c) the condensate fraction /o, and (d) the oscillator 
length parameter a for nonrelativistic boson systems with iV = 2000 and TV = 250 in a static equilibrium with a r rms = 10 fm, 
plotted as a function of temperature. 



What is the depth of the mean-field potential that holds the pions together in static equilibrium for a given r rms ? 
The depth of the potential is approximately huj(r rms /a) 2 /2 (see Eq. (1)). For r mrs = 10 fm and T ~ 120 MeV, the 
results in Figs. 8(a) and 8(c) show that the depth of the potential needs to be about 18(MeV) x 2.5 2 /2 ~ 56 MeV 
for N = 2000, and is about 25(MeV) x 3 2 /2 ~ 112 MeV for N = 250. These are not very deep potentials. It will be 
of interest to determine theoretically the mean- field potential for an assembly of pions at different temperatures. 

We reach the following conclusion from the above study: if a nonrelativistic pion system maintains a static equilib- 
rium within its mean field, and if it contains a root-mean-squared radius, a pion number, and a temperature typical 
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of those in high-energy heavy-ion collisions, then it will contain a large fraction of the Bose-Einstein pion condensate. 
The pion condensation will affect the parameter A in momentum correlation measurements. 

The evolution of pions in high-energy heavy-ion collisions involves dynamical motion and may not be in a state of 
static equilibrium. The static solutions examined here serve as supplementary tools relative to which the effects of 
the dynamical motion and non-equilibrium effects may be investigated. 



IX. RELATIVISTIC TREATMENT OF A BOSON GAS IN A HARMONIC OSCILLATOR 

For pions in the environment of a high-energy heavy-ion collision, the pion temperature is of the order of the pion 
rest mass and a relativistic treatment of the pion motion is needed. We therefore examine a boson in an external field 
characterized by a time-like vector interaction A (r), a space-like interaction A(r) and a scalar interaction S(r). The 
Klein-Gordon equation for the motion of the boson is 

[p - A (r)} 2 - [p - A(r)] 2 - [m + S(r)Au(r) = 0. (62) 



Different types of interaction potentials will lead to different single-particle spectra and different Bose-Einstein con- 
densations that will need to be explored in more detail in the future. We shall examine here at this stage only 
the simplest kind of exactly solvable potential that is closely connected to the harmonic oscillator potential in the 
nonrelativistic limit. Accordingly, we study scalar interactions S(r) and introduce the interaction interaction V(r) 
related to S(r) by 

V(r) = S(r) + M! (63) 

The V(r) and the S(r) potentials approach each other in the nonrelativistic limit of m — > oo. In terms of V(r), we 
have 

[m + S(r)] 2 =m 2 + 2my(r), (64) 
and the eigenvalue equation for relativistic motion with only a scalar interaction becomes 

+ v(r) } u(r) = ^tr" u(r) s eu(r) ' (65) 

where the eigenvalue e is related to the particle energy p n by 

Pa = E= Vm 2 + 2me. (66) 

To make the problem simple and to connect with earlier exactly solvable nonrelativistic solutions, we choose to 
consider V(r) to be the same harmonic oscillator potential of Eq. (1), 

V{v) = imwV. (67) 
The eigenenergy of the relativistic boson is exactly soluble and is 



E n = yjm" 1 + 2me n , (68) 

where 

e n = {n+^)hj. (69) 
We likewise introduce the recalibrated E n measured relative to the energy of the n = state 



E n = y to 2 + 2m(n + -hw) - y m 2 + 2m x -fku. (70) 

Instead of the nonrelativistic condition of Eq. (6), the relativistic condensate configuration condition becomes 

°° —pit 

N = N + N T = ^- + Y. ? nZe Te ■ W 
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FIG. 9: (Color online) The fugacity parameter z satisfying the relativistic condensate configuration condition Eq. (71) for 
different boson numbers N and huj/m — 20/140, as a function of the temperature T/fko. Fig. 9(b) is an expanded view of Fig. 
9(a) in the z ~ 1 region. 

where g n — (n + l)(n + 2)/2. 

To solve for the fugacity parameter z in the relativistic case, it is necessary to specify fko/m, the ratio of harmonic 
oscillator energy scale fko to the rest mass m of the boson. A small fko/m ratio approaching zero corresponds to the 
nonrelativistic limit and a large ratio relative to zero leads to the relativistic case. 

We are interested in the case where the mass m of the boson is of the order of the gas temperature T. We shall 
see in Fig. 13a below that a boson system with N = 2000 pions, T ~ m, and a root-mean-squared radius of 10 fm 
corresponds to a harmonic oscillator energy fko as 40 MeV which is a substantial fraction of the rest mass m. We shall 
therefore investigate relativistic boson systems with fko/m = 40/140 in our numerical studies. 

With this specification of fko/m while N and (3fko = fko/T are held fixed, the relativistic condensate configuration 
condition (71) can be solved numerically to determine the unknown z. We show in Fig. 9 the fugacity z which satisfies 
the relativistic condensate configuration condition for different temperatures T/fko and different boson numbers N. 
To get a better view of the z values, we show an expanded view of Fig. 9(a) in the z ~ 1 region in Fig. 9(6). 

We observe that the fugacity parameter z is close to unity in the strongly coherent region at low temperatures. 
Upon a comparison of Fig. 9 with Fig. 1, one notices that the shapes of 2 as a function of T/fko for the relativistic 
and the nonrelativistic cases are very similar, except that the scale of the temperatures are much reduced for the 
relativistic case. For N — 2000, the condensate temperature occurs at T/fko ~ 3.5 in the relativistic case, in contrast 
to the nonrelativistic case at T/fko ~ 11. For N = 250, the condensate temperature occurs at T/fko ~ 1.9 in the 
relativistic case, in contrast to the nonrelativistic case of T/fko ~ 5.1. To see why these large changes occur, we note 
that the recalibrated energy expanded in powers of 1/m is 

^_ [(n + 3/2) M ' + ^ (72) 
Zm m z 

For the relativistic harmonic oscillator potential we have chosen, the spectrum of E n is nearly the same as those in 
the nonrelativistic case of e n for small values of n. However, the spectrum for large values of n is greatly compressed 
by the presence of the second term with a negative sign in the above equation. As a result, a large number of chaotic 
particles can be accommodated even at a lower temperature, leading to a large shift of the condensate temperature 
in units of fko in Fig. 9 when relativistic effects are included. 

After the value of the solution z is obtained, No = z/(l — z) and Nt can be subsequently determined to give the 
condensate configuration specified by the condensate fraction /o and the chaotic fraction fa. 

The condensate fractions fo(T) as a function of T/fko calculated with the fugacity parameters of Fig. 9 for dif- 
ferent boson numbers N and fko/m = 40/140 are shown as different curves in Fig. 10. The abscissa labels for the 
corresponding chaotic fraction fa{T) = [1 — fo(T)] are indicated on the right. We observe that the behavior of the 
condensate fraction in the relativistic case is similar to the nonrelativistic case, with the exception of the shift of 
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FIG. 10: (Color online) Different curves represent condensate fractions fo(T) as a function of T/hui for different boson numbers 
N, calculated with the relativistic condensate configuration condition Eq. (71) for hu)/m — 40/140. The abscissa labels for the 
corresponding chaotic fraction fr(T) = [1 — /o(r)] are indicated on the right. 

temperature T/Hlu to lower values. Again, the transition from the condensate phase to the chaotic phase occurs over 
a large range of temperatures and is therefore not a sharp first-order type transition. The complementary chaotic 
fraction fr{T) increases gradually as the temperature increases, reaching the value of unity at large T/hcv. 

X. SPATIAL AND MOMENTUM DISTRIBUTIONS IN THE RELATIVISTIC CASE 




r/a , pa/h 



FIG. 11: (Color online) The spatial density distribution p r (r) in units of a -3 , expressed as a function of r/a, and the momentum 
density distribution p p (p) in units of (a/fi) 3 , expressed as a function of pa/h, for the relativistic case of hw/m = 40/140 with 
N = 2000 at different temperatures. 

The knowledge of the fugacity parameter for different temperatures allows one to determine the occupation numbers 
at different single-particle states. These occupation numbers and the absolute square of the single-particles wave 
functions give the spatial and momentum densities of the system at different temperatures. As we remarked previously, 
p p (p) and p r (r) have the same shape when properly scaled. It suffices to consider the spatial density p r (r) given by 

oo -/3E n 

Pr (r) = GW(r,r) = J2 - <(rK(r), (73) 

n— 
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where n represents the set of quantum numbers {n r lm} of a harmonic oscillator state, u n (r) is the harmonic oscillator 
wave function normalized to J dr\u n (r)\ 2 = 1, 

u n (r) = N nrl x l e~ x2 l 2 L l +hx 2 )Y lm {6^), (74) 

n = 2n r +1, x = r/a, L nr 2 (x 2 ) is the associated Laguerre polynomial, and 




TV ftro 

FIG. 12: (Color online) The root-mean-squared radius in unit of a and the root-mean-squared momentum in units of ft/a, as 
a function of T/hcj for different numbers of bosons in the system. 

We plot in Fig. 11 the spatial and momentum distributions of the system with N = 2000 for the relativistic case 
of fko/m = 40/140 as a function of their dimensionless variables r/a and pa/h respectively. One observes that up 
to T/fku ~ 3 the system has a small spatial or momentum size and there is a substantial condensate fraction in the 
system. In Fig. 12 we plot the root-mean-squared radius in unit of a, r lms /a = a/ ((r/a) 2 ), and the root-mean- 
squared momentum in unit of H/a, p ruls a/h — \J ((pa/h) 2 ) , as a function of T/huj. For N = 2000, the quantity 
frms/a is slightly greater than 1 up to T/hui ~ 2, and it increases relatively rapidly to about 6.5 at the condensate 
temperature, T c /hu> ~ 3.5. It increases at a relatively slower rate at temperatures above T c . 

We can carry out an analysis to inquire the following: If a system of N = 2000 relativistic pions is held together 
by its mean field, taken to be a harmonic oscillator, and if it comes to a state of static equilibrium with a given root- 
mean-squared radius r rms at a temperature T, what is the condensate fraction of such a system? We shall therefore 
examine a pion system with r rms = 10 fm, a temperature range from 80 to 160 MeV, and the number of pions to 
be N = 2000 (for a central RHIC Au-Au collisions at ^/sjvjv = 200 GeV). We have chosen frw/m = 40/140, to be 
approximately self-consistent with the value of fkv extracted form such an analysis (see Fig. 13(a)). 

If the value of r rms is fixed as given, the quantities hio and T are then related by the set of parametric equations 
of (60) and (61). Using the function r lms /a = F^(x) of Fig. 12 in the above equations and varying x = T/hw for a 
fixed r rms , one can determine the energy hui as a function of T for the case of N = 2000. The results are shown in 
Fig. 13(a). One finds that for the pion system with a given root-mean-squared radius of 10 fm, the value of hu> ranges 
from about 30 to 53 MeV for N = 2000, with an average of about 42 MeV. The ratio of T/fiuj is about 2.7 to 3 for 
N = 2000 as shown in Fig. 8(6). From these ratios of T/hcu, one can use Fig. 2 to find out the condensate fraction. 
The condensate fractions fo(T) for a pion gas at various temperatures with N = 2000 are shown in Fig. 8(c). One 
finds that fo(T) is between 0.5 and 0.7 for N — 2000. The knowledge of Hui in Fig. 8(a) allows one to determine the 
values of a as a function of temperature as shown in Fig. 8(d). The oscillator length a is between 3 and 2.2 fm for 
N = 2000. 

The relativistic analysis indicates that the relativistic effects change the single-particle spectrum and shift the 
locations of the condensate fraction in units of fko. The condition of maintaining a system size with a root-mean- 
squared radius of 10 fm recalibrates and raises the oscillator energy hu for the relativistic case, as compared to the 
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FIG. 13: (a) the potential strength hu, (b) the ratio T/hai, (c) the condensate fraction fo, and (d) the oscillator length parameter 
a for relativistic boson systems with N = 2000 and N = 250 in a static equilibrium with a r rms = 10 fm, plotted as a function 
of temperature. 



nonrelativistic case. As a consequence, the condensate fraction for N = 2000 is modified from / ~ 0.67 to 0.8 in the 
nonrelativistic to fo ~ 0.5 to 0.7 in the relativistic case. There is a small reduction of the condensate fraction, but 
the condensate fraction remains quite large in the relativistic case. 

We again reach the following conclusion from this study: If a relativistic pion system maintains a static equilibrium 
within its mean field, and if it contains a root-mcan-squarcd radius, a pion number, and a temperature typical of those 
in high-energy heavy-ion collisions at RHIC, then it will contain a large fraction of the Bose-Einstein pion condensate. 



XI. DISCUSSIONS AND CONCLUSIONS 



As the chaoticity parameter A has been widely used in all HBT measurements, we are therefore motivated to 
investigate an exactly solvable problem to study the momentum correlation function for a noninteracting boson gas 
assembly held together in a harmonic oscillator potential at various temperatures. In the process, we find that 
the phase transition from the Boson-Einstein condensate to the chaotic phase occurs gradually over a large range 
of temperatures, with the condensate fraction fo(T) varying approximately as 1 — (T/T c ) 3 , where the condensate 
temperature T c is approximately given by (7V/1.202) 1 / 3 fiij. The spatial and the momentum radii of the system are 
small in a condensate at low temperatures, of the order of a few oscillator units, increasing in size as the temperature 
reaches the chaoticity limit. 

From the momentum correlation function, we can determine the \(p,T) parameter and the HBT radius -Rhbt- 
We find that the X(p,T) parameter is a sensitive function of both the pair momentum p and temperature T. For 
a temperature above the condensate temperature, the X(p, T) parameter is 1 for all momentum p. However, for 
temperatures below and even substantially below the condensate temperature, X(p, T) is small and close to zero for 
small pair momentum p, but it increases and saturates at X(p,T) = 1 at large pair momentum p. The location 
where X(p, T) attains unity changes with temperature. The lower the temperature, the greater is the value of p at 
which the X(p, T) attains the value of unity. Because the X(p, T) parameter attains the value of unity for systems at 
temperatures much below the condensate temperature, the occurrence of A = 1 is not consistently correlated with the 
absence of a condensate fraction. Only in the region of small p will the parameter X(p, T) be correlated with, but not 
equal to, the chaotic fraction fr(T) of the system. 

We find that the HBT radius -Rhbt increases gradually with increasing pair momentum p and temperature T. 
However, for small value of p, the HBT radius increases only slowly with increasing temperature at low tempera- 
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tures and it then increases rapidly and abruptly as the temperature approaches the condensate temperature. The 
temperature dependence of the HBT radius at small p values correlates well with the temperature dependence of the 
root-mean-squared radius of the system. 

It is of interest to inquire the degree of coherence of pion systems produced in high-energy heavy-ion collisions. 
We have examined both cases of pions as a nonrelativistic and a relativistic gas in a harmonic oscillator potential. 
If a pion system maintains a static equilibrium in its mean field, and if it contains pion numbers from N = 250 to 
N = 2000, a temperature in the range from 80 to 160 MeV, and a root-mean-squared radius of 10 fm, (typical of 
those one encounters in high-energy heavy-ion collisions), then it will contain a large fraction of the Bose-Einstein 
pion condensate. While the details of the nonrelativistic and relativistic calculations are presented in Sections VIII 
and X, we can provide simple arguments here to indicate that these are reasonable results based on plausible physical 
principles. Bose-Einstein condensation occurs when the temperature is below the condensate temperature, which is 
a few units of two. We need to estimate the energy scale fkj for the pion system. The energy scale can be estimated 
by knowing the length unit a. One expects that the pion system with a root-mean-squared radius of 10 fm would be 
contained within a few units of this length a, leading to a rough estimate of the length unit a to be about a few fm. 
By dimensional analysis, the energy scale associated with this length unit a for a pion is fiw = tr/ma , which gives a 
value many tens of MeV for two. With a temperature of T = 120 MeV or 140 MeV, we obtain the ratio T/tuv of a few 
units which would correspond to a T/fkv ratios with a substantial condensate fraction. We can therefore understand 
that the occurrence of the pion condensation in static equilibrium arises because the pions are massive particles, and a 
large number of pions are produced and concentrated in a small spatial volume characterized by a root-mean-squared 
radius of only 10 fm. The pion gas in static equilibrium is therefore in the realm of low-temperature boson systems 
with possible occurrence of Bose-Einstein condensation. 

The evolution of pions in high-energy heavy-ion collisions contain dynamical motion and may not be in a state 
of static equilibrium. How the dynamical motion of the pions will modify the coherence of the system will be an 
interesting subject for future investigations. 

While we await future theoretical investigations, it is of interest in the meantime to discuss possible modifications 
of the static results obtained here in the presence of a collective expansion. One expects that the collective expansion 
will not alter the energy ordering of the states of the system. Because the Bose-Einstein condensation depends on 
the relative ordering of the energy of the states, the coherence may not be greatly affected. However, the average 
pair momentum in the direction of the expansion will be greatly boosted. The X(p, T) parameter for the expanding 
coherent source as a function of the pair momentum would likely retain a shape similar to the static case, but with 
the p boosted by the collective expansion. In this connection, it is interesting to note that the experimental A values 
plotted as a function of the pair transverse momentum has a shape [54, 55] quite similar to the shape of the X(p, T) 
parameter plotted as a function of the average pair momentum pa/H in Fig. 6(a). Although alternative explanations 
in terms of a decrease in the resonance decay contributions at higher px have been presented, it will be of interest to 
explore whether the behavior of A as a function of the pair transverse momentum may be due to the occurrence of 
an expanding Bose-Einstein condensate. 

The theoretical HBT radius for a static source increases slightly as the pair momentum increases, while the exper- 
imental measurements gives an HBT radius decreasing as px increases. The theoretical HBT radius may be more 
sensitively affected by the expansion dynamics because the collective expansion boosts not only the average pair 
momentum but also the relative momentum between the correlated pair, the boost being the greater the larger the 
magnitude of the pion momentum. As a larger relative momentum leads to a smaller HBT radius, the HBT radius 
therefore decreases as a function of the pair momentum. Clearly, whether future analyses bear out this possibility 
will be of great interest. How the collective pion motion will affect quantitatively the HBT radius for a boson system 
with varying degrees of coherence is therefore an interesting subject for future investigations. 

It has been proposed that the question of whether an observation of A < 1 is due to coherence or due to contamination 
from particles from far outside the source volume can be tested by analyzing three-particle correlations [56]. Such 
analyses of data at both SPS and RHIC have been consistent with the chaotic conjecture [57]. However, as we note 
that the A(p, T) parameter can assume the value of unity in certain kinematic regions even for significantly coherent 
systems with a temperature much below the condensate temperature, the attainment of A = 1 cannot be a unique 
signature of the chaoticity of a system. However, how the coherence of the boson system may affect three-body 
correlations has not been worked out explicitly and merits further investigations to clarify the situation. 

It needs to be emphasized that to make the problem tractable as an exactly solvable model, we have specialized to 
a static treatment of the boson system in equilibrium in both a nonrelativistic and a relativistic harmonic oscillator 
potential, for which analytical eigen-energies and eigenfunctions can be readily available. This is a simple model of 
noninteracting boson gas in an external potential without two-body interactions. Even with such an idealization, a 
wealth of new information on the coherence and two-particle momentum correlation functions as well as the chaoticity 
parameter has been obtained as a function of the attributes of the boson environment. 

The harmonic oscillator potential introduced here can arise from an external trap, as in atomic physics. In high- 
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energy heavy-ion collisions, the harmonic oscillator potential can arise approximately from the mean-field potential 
experienced by a pion, due to the interactions generated by other pions and medium particles. While approximating 
the pion mean-field potential as a harmonic oscillator potential can yield gross features and a wealth of information, 
a more accurate determination of the pion momentum and correlation functions will require a better description of 
the pion mean- field potential. As the mean- field potential depends on the pion density as in the Glauber theory [53], 
and the equilibrium pion density depends in turn on the mean-field potential, it will be necessary in the future to 
study the pion mean-field potential and the density self-consistcntly in a pion condensate. Besides these mean-field 
interactions between pions, the remaining residue interactions will give rise to additional complications which may be 
studied in the future. 

We have described the correlations in a static equilibrium environment and we need to discuss how the time 
dependence can be handled. An accurate dynamical treatment will examine the time evolution of the system for a set 
of given initial conditions, as in a time-dependent Hartree approximation following the technique of the time-dependent 
Hartree-Fock approximation developed in nuclear physics [58, 59] 

While we await future work on the time dependence of the correlation function, the static results obtained here can 
be used for experimental comparison if the time dependence of the external field is such that the external potential 
is suddenly removed, as in a typical condensed matter experiment with trapped atomic particles in a condensate. 
In these low-temperature measurements with atoms, the trapped atoms before being released are now described as 
having an equilibrium momentum distribution in momentum space, appropriate for the system in a given external 
field at a given temperature. The sudden removal the external field allows the initial momentum distribution of 
the particle to be frozen at the moment of the external field removal, as appropriate under the application of the 
sudden approximation in quantum mechanics. Subsequent free streaming of the particles without the external field 
and mutual interactions (except for the additional correction of the gravitational field or other extra forces applied to 
the particles) allows the reconstruction of the momentum distribution of the source at the moment of its freezing out. 
In measuring the arrival times and arrival positions of the particles of a correlated pair in Ref. [50] , the quantities 
that are in effect measured are the momenta of correlated pairs from which the average momenta and the relative 
momenta of the pair can be collected and examined. The perspectives of studying the correlation in momentum space 
presented here offer useful complementary viewpoints to the theoretical and experimental works in atomic physics 
have been focused so far on the correlation function in configuration space. 

The results obtained here can be approximately applied to heavy-on collisions if the explosive expansion is so rapid 
that it can be approximately described as a sudden removal of the external field. In that case, the static initial 
momentum distribution and correlations of the particles would be frozen at the moment of the external field removal 
and show up as particles reaching the detectors by free streaming. In this respect, it is of great interest to examine 
in the future a dynamical model of the expansion of the pion gas and study how the explosive expansion will affect 
the momentum correlation function. 
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